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Abstract 



In paper I [M. Znojil and G. Levai, Phys. Lett. A 271 (2000) 327] we intro- 
duced the Coulomb - Kratzer bound-state problem in its cryptohermitian, 
T'T— symmetric version. An instability of the original model is revealed here. 
A necessary stabilization is achieved, for almost all couplings, by an unusual, 
negative choice of the bare mass in Schrodiner equation. 
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1 Introduction 



Intuitively one feels that for Schrodinger equations 
r L{L + 1) 



2m 



2 



^(x) + V{x) ^{x) = E ^{x) (1) 



there should exist a close connection between the reality of potential V{x) 
and the reality of the corresponding energies E. Unfortunately, such a type 
of intuition proves deceptive. Recent studies (e.g., or ^) showed that 
many manifestly non-Hermitian potentials, e.g., 

^(ss)(a;) =x2(ia;)^^ 6>0 (2) 

still lead to a full reality of the spectrum. The key to such an unexpected 
phenomenon can be seen in the Bender's and Boettcher's [T] fortunate choice 
of a complex integration contour x = x^^^\s) in eq. (JT]). Its asymptotes 

a;(«^)(.) ^ 1 (3) 
I se-"^ , s -C -1 

were restricted to the 5— dependent interval of angles, 

n f n 3n \ 
+ o e — — ^, — — ^ . (4) 



2 \^4 + 45'4 + 45 

The curve itself was required not to cross the singularity of V^^^\x) in the 
origin, ix''^^''(0) > 0. In this setting one can impose the standard Dirichlet 
boundary conditions at the ends of the left-right-symmetric curve of complex 
coordinates, \l/(x^'^^^(±oo)) = 0, with the computationally preferred slope 
lying precisely in the center of the interval, 

^ 2 + 26 ~ 2 ■ ^ ^ 

In ref. [T] it has been emphasized that potentials ([2]) as well as paths of x and 
angles ([5]) were chosen as symmetric with respect to the combination of the 
parity-reversal symmetry mediated by the operator V with the time-reversal 
symmetry represented by operator T (cf. also ref. [3J in this respect). In 
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ref. |1] it has been added that for the other ehgible domains of angles, say, 
for 

^ ^ Stt 57r \ 

(6) 



TT 



57r 



2 \^4 + 4(5 4 + 4(5^ 
the reahty of the spectrum breaks down at some non-vanishing exponents 
6 < 6q. In this sense the specific PT— symmetric choice of ([2]) - dl]) giving 
Sq = may be considered optimal. 

The discussions in refs. [H H] did not involve the negative exponents 
S and, in particular, the short-range models where V{oo) = 0. The gap 
has partially been filled by ref. [5] where we studied eq. ([I]) with one of the 
simplest possible asymptotically vanishing PT— symmetric potentials of the 
Coulomb-Kratzer two-parametric form. 



V^(x) = y(^^)(x) = - + 4- 



X 



(7) 



This model admits ip G (0, vr) (cf. eq. 



with 26 = —1). From eq. 



giving (p^^^^ = tt/2 one arrives at the U-shaped complex-coordinate contours 
r^{BB) ga]-Qpig(j in Figured] where the cut is assumed from a; = upwards. 
Marginally let us emphasize that our Schrodinger eq. ([1]) in the most common 
physical setting using an integer angular momentum £ = 0, 1, . . . should in 
fact be considered with the "centrifugal-like" term of a non-integer strength 
L{L + l) = i{i+l) + F in general. 




Figure 1: The optimal, U-shaped contour of the complexified coordinates 
x^^^\s) for the Coulomb-Kratzer PT— symmetric potential ([7]). 
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At the time of the pubhcation of ref. [5] (to be cited as paper I from now 
on) the physical meaning of the similar models remained still rather obscure. 
Many authors studied and interpreted them as mere effective non-Hermitian 
simulations of spectra, not allowing any immediate physical interpretation 
of the related wave functions \E'(x) G 1L2{M). Although we also accepted 
the same philosophy in paper I, we were aware of the fact that such an atti- 
tude significantly weakened the impact and practical applicability of similar 
studies. 

Fortunately, the subsequent development of the subject clarified that the 
potentials as exemplified by eq. ([7j) can be interpreted as fully compatible 
with the standard postulates and probabilistic interpretation of Quantum 
Mechanics. One of the most straightforward mathematical keys to the reso- 
lution of such an apparent puzzle can be seen in the existence of a suitable 
non-unitary inveitihle map Q between some manifestly non-Hermitian Hamil- 
tonians H ^ W and their manifestly Hermitian partners f) = VLHQt^ (cf., 
e.g., ref. [6j| for a compact explanation of this mathematical idea). 

From the point of view of physics, the historical origin of the idea of 
relevance of isospectrality between f) and H can be traced back to the study of 
models of atomic nuclei [7j. There an explicit example of operator Vt ^ (^^)^^ 
has been provided by the generalized Dyson mappings [8] . Our recent return 
to these physical studies in our mathematical review [6] showed that for the 
"PT— symmetric models all the probabilistic physical postulates of quantum 
theory remain valid. 

Among immediate and most recent phenomenological applications of non- 
Hermitian, "PT— symmetric operators H ^ with real spectra let us men- 
tion here just the preprint [9] dealing with a PT— symmetric version of a flat 
Friedmann model in quantum cosmology. In such a broader physical context 
we feel particularly inspired here by one of technical questions discussed in 
this paper and concerning the possible instabilities of generic PT— symmetric 
systems. In this sense we also returned to our older results of paper 1 which 
will be re-evaluated, corrected and re-interpreted in what follows. 
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2 Free motion along asymptotes 



In the majority of presentations of Schrodinger eq. ([T]) in textbooks one works 
with the real x specifying the position of a particle or quasiparticle which 
carries a constant mass m = mo > 0. The influence of external forces is 
modeled solely by a potential V{x). During the last few years a manifest 
coordinate-dependence of the mass term has been allowed as well [10]. The 
choice of m = m(x) opened new perspectives in an optimal description of 
the effects of medium. 

This idea could easily be transferred to the present class of models where 
X = x{s) is complex and where the effect of the potential becomes negligible in 
the asymptotic domain of |s| ^ 1. In such a setting the mass can be perceived 
as a potentially position-dependent complex quantity, m = m[x{s)] G (D. 

At the large |s| our Hamiltonians get approximated by the kinetic-energy 
operator T which, by itself, gets complexified in the light of eq. ([3]), 



Once we introduce the asymptotically constant complex effective local mass 
meff[x{s)] it will only depend on the slope ip and on the sign of s, 



This observation is too abstract, for several reasons. First of all, a sub- 
tle balance between the left and right branches of wave functions \I'[x(s)] 
exists and reestablishes the reality of the energies for numerous complex in- 
teractions V"[a;(s)] [IT]. Secondly, for m = m(x) the well-known von Roos' 
jlOj ambiguity of the kinetic energy would emerge at the finite values of s. 
For complex x{s) the manifest introduction of the coordinate-dependence in 
the mass might also lead to many other technical complications. For these 
reasons our present attention will solely be paid to the models where me// 



T 



2mo dx"^ 





T 




(9) 
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remains constant. Using just the asymptotically vanishing potentials exem- 
plified by eq. ([7]) and assuming the local reality of the kinetic energy we shall 
only make a choice between cp = and = 7r/2. In this way we encounter 
either the entirely traditional textbook straight-line models at (/9 = or their 
U-shaped-line innovations at = 7r/2. 



Im X 


x(s) 




J Re X 



Figure 2: The e = contour x^'^\s) = xL (s). 



As long as the former case is very traditional let us only discuss the choice 
of ip = 7r/2 giving the U-shaped contours sampled in Figure [H Since both 
their asymptotes parallel the upper imaginary half-axis (i.e., a cut from x = 
upwards), the phase of the complex numbers will be assumed lying in the 
interval (— 37r/2, 7r/2). Under such a convention and in terms of a suitable 
width parameter e > we may parametrize the contours of Fig. [1] as follows, 

{-i(s + le) - £, s G (-00, -|£:), 
^gi{s/.-i/2.)^ 3e(-f£,f£), (10) 

In the complex plane of x the latter curve exhibits the double-refiection left- 
right symmetry x{—s) = —x*{s) which combines the spatial refiection V 
with the complex conjugation T (let us recollect that T : i — > — i mimics 
time-reversal). Our next Figure [2] shows how such a curve of the complex 
coordinates could be deformed in the limit 6 = 0. It still encircles the origin 
at a distance but its asymptotes already strictly coincide with the upper 
imaginary half-axis. 
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Let us emphasize that for Lp = 7^/2 the coordinate- independence of the 
effective mass simphfies the kinetic-energy operator 

Surprisingly enough, it acquires the wrong sign in the sense that its spectrum 
becomes unbounded from below at the positive "bare mass" mo > 0. This 
would make the whole system unstable with respect to small perturbations 
and, hence, useless for any phenomenological purposes. 

There are hints that also in a field theoretical framework similar consid- 
erations hold concerning negative kinetic energy encountered during quan- 
tization of classical phantom Lagrangians [9]. This encourages us to make 
our argument more quantitative. Let us recollect the asymptotic form of our 
Coulomb - Kratzer Schrodinger equation at |s| ^ 1, 

-^,n^)=a-En^)^ a^ = ^>0. (12) 

Distinguishing between the positive-energy domain {E = k"^ > k > Q) 
and the negative-energy domain [E = —k^ < 0, k > 0) we may employ the 
general superposition formula 

^(x) = C+ + C_ , X = x{s) ~ i |s| + . . . 

and insert the pair of linearly independent solutions of eq. (fT2!) . 

giiafcx _ ^Tak\s\ ^ E = k"^ > , 



^±{x) OC 



e 



±aKX _ g±ia«|.| E = -K^ >0. 



The upper option proves linked to the asymptotically vanishing bound states 
which were constructed in paper I at > 0. In parallel, the lower-line option 
reveals the admissibility of the free plane-wave states at all the negative 
energies. This implies, in a way unnoticed in paper I, that there exists also 
a continuous part of the spectrum which remains unbounded from below. 

In the light of the latter semi-intuitive argument our bound-state model 
of paper I appears unstable with respect to perturbations and, hence, deeply 
unphysical. This forced us to write the present addendum to paper I showing, 
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in essence, that a complete remedy of such a very serious shortcoming is 
unexpectedly easy. Our key idea is that once we deform the coordinates we 
must also turn attention to the underlying theory (cf. [6j) and re-analyze all 
the questions of the mathematical consistency of the model. 

3 Amended Coulomb - Kratzer bound states 

First of all, we must impose the forgotten but essential requirement of sta- 
bility, i.e., of the boundedness of the spectrum from below. In this sense 
our present main result is that the latter requirement can be satisfied rather 
easily. Formally, it appears equivalent to the reversal of the sign of the bare- 
mass parameter, mo = —mi < 0. In order to explain this usual amendment 
of the model let us first replace eq. f|T2l) by the modified asymptotic equation 

-^,^{x) = -P^E^{x), /?^ = ^>0 (13) 

which, mutatis mutandis, implies that 

g±/3fcx _ g±i/3fc|s| ^ E = k^ >0, 

Using the same argument as above we deduce that the continuous spectrum 
is positive and that the discrete bound-state energy levels may be expected 
negative. In this way the structure of the spectrum of the non-Hermitian 
Coulomb-Kratzer model of paper I is thoroughly modified and made more 
similar to its well known textbook Hermitian-Coulomb-Kratzer predecessor. 

We saw that the spectrum of our particular illustrative example as well 
as of all the similar asymptotically non-interacting models may be made 
acceptable, on physical grounds, only if we complement the complexification 
of coordinates by the parallel adaptation of the bare mass. We must keep in 
mind that even the complexification of x{s) itself is often perceived as unusual 
since it causes the complete loss of the observability of coordinates. This step 
has only recently been accepted as an admissible innovative model-building 
recipe which characterizes almost all PT— symmetric quantum models. 
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Our present key recommendation of the choice of a negative mass tuq may 
look equally counterintuitive. We believe that it deserves to be accepted 
on similar grounds, as a mere very natural mathematical consequences of 
the complexification of x{s). Indeed, the complexification of x{s) immedi- 
ately implies a breakdown of the traditional split of the Hamiltonian into 
its kinetic- and potential-energy parts so that the switch to the negative 
value of the bare mass parameter m = mo = —mi < is in a one-to-one 
correspondence with the guarantee of the stability of the system in question. 

Let us return to a constructive demonstration of consistence of the negative- 
mass bound-state problem, recollecting first the results of paper I where 
the solvable Coulomb-Kratzer potential has been inserted in the traditional, 
positive-mass Schrodinger equation. 



L{L + 1) .Z 



^(t) = ^^(t). (14) 



This equation has only been considered at non-integer L in paper I. Here, 
we shall accept the same constraint and assume that L ^ 0, ±1, . . .. This 
enables us to retype the formula for the discrete eigenvalues from paper I, 

2 



E. 



±n 



n = 0,l,.... (15) 



_2L + l±(2n + l)_ 

One feels puzzled when seeing that all of these eigenvalues are positive. 
Indeed, the negative bound-state energies would be generated by the real 
Coulomb and Coulomb - Kratzer potentials [12] . 

In the light of our preceding considerations we know that the spectrum 
f|T5|) must be discarded as unstable. This resolves the latter paradox and, 
marginally, it also could throw new light on some recent attempts of using 
the Coulomb-like complexified potentials and/or the negative-mass option 
in different contexts [I3l [lU |T5]. For example, Mostafazadeh ^16\ noticed 
that in the latter preprint [15] the PT— symmetry violation caused by the 
complex-scaling transformation of x has led to a negative bare mass but 
disabled the authors to cope properly with boundary conditions. Actually, 
the similarity transformation expressing the complex scaling transformation 
of X violates the VT symmetry and either introduces an imaginary part in 
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the energy or transforms normalizable states into non-normalizable ones and 
vice versa |T7j. 

The mathematical core of our present proposal is different. In essence, 
our present recipe degenerates to the mere change of the overall sign of the 
tentative Coulomb Hamiltonian of ref. [5]. The corrected negative- mass ver- 
sion of the present update of the PT— symmetric Schrodinger equation ([1]) 
in its Coulomb - Kratzer exemplification becomes obtainable from eq. (fT4l) 
by formal substitution E —E. Although this implies the reversal of the 
sign of Z, such a modification of the potential is inessential since the eigen- 
values f|T5l) themselves are only proportional to Z^. Our final negative- mass 
Coulomb - Kratzer Schrodinger equation may be written in the form 

L{L + 1) 



Z 
I — 

X 



X G x)J-^ \s) 



(16) 



yielding the bound-state-energy formula 

Z 



E = E. 



±n 



H 2 



n 



0,1, 



(17) 



_2L + 1 ± {2n + l] 

In Figure [3] this coupling-dependence of the energy levels is illustrated via 
the ten lowest bound states. 
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Figure 3: Spectrum of —Kn,a at Z = e = 1 as a function of 2L -|- 1 with 



n = 0, 1, . . . , 4 and a = ±1. 



Let us summarize that after we changed the sign of the bare mass the spec- 
trum of our amended PT— symmetric Coulomb-Kratzer interaction model 
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looks qualitatively similar to its Coulomb and Coulomb-Kratzer Hermitian 
predecessors. Its continuous part of the spectrum is "well-behaved" and non- 
negative, i.e., it is bounded from below - this guarantees the stability of the 
system. Similarly, all the discrete energy levels only posses single accumula- 
tion point at -E = 0. 

Still, the differences illustrated by Figure [3] are also worth mentioning 
(cf., e.g., [18] for comparison). First of all, in contrast to the Hermitian 
Coulomb-Kratzer model the present discrete spectrum is composed of the 
two qualitatively different families of levels which are distinguished by the 
±— ambiguity in formula f|T71) . As a consequence, the traditional "fall of 
the particle on the center" known from the textbooks [T2| [18] is now re- 
peated at any integer "singular value" of our Kratzer-coupling-dependent 
non-integer parameters L = L{F) = 0, 1, . . .. This observation also offers a 
purely physical explanation why we had to omit these singular values from 
our considerations. 

In place of the picture we may also employ the following reparametrization 
oi 2L + 1 = Mq + cos^ a > where the integer part Mq > of this parameter 
is complemented by a small positive residuum cos^ a < 1 where a G (0, 7r/2). 
This decomposition of L = L{Mo, a) leads to the compactification of the 
ground-state-energy formula 

1i 

mm(sm a, cos^ a) 

Although this function of a represents the lower bound of the whole spec- 
trum, this function is, by itself, unbounded from below. This means that 
in the "allowed" vicinity of the "excluded" limiting values of a = and 
a = tt/2 our system still gets very strongly bound. Moreover, even quite 
far from a = and a = tt/2 all the low-lying spectrum remains extremely 
sensitive to the small perturbations or variations of the coupling constant F. 
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